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ABSTRACT. We consider a general family of regularized Navier-Stokes and Magneto- 
hydrodynamics (MHD) models on n-dimensional smooth compact Riemannian mani- 
folds with or without boundary, with n > 2. This family captures most of the specific 
regularized models that have been proposed and analyzed in the literature, including 
the Navier-Stokes equations, the Navier-Stokes-a model, the Leray-a model, the Mod- 
ified Leray-a model, the Simplified Bardina model, the Navier-Stokes-Voight model, 
the Navier-Stokes-a-like models, and certain MHD models, in addition to representing 
a larger 3-parameter family of models not previously analyzed. This family of models 
has become particularly important in the development of mathematical and computa- 
tional models of turbulence. We give a unified analysis of the entire three-parameter 
family of models using only abstract mapping properties of the principal dissipation and 
smoothing operators, and then use assumptions about the specific form of the parameter- 
izations, leading to specific models, only when necessary to obtain the sharpest results. 
We first establish existence and regularity results, and under appropriate assumptions 
show uniqueness and stability. We then establish some results for singular perturba- 
tions, which as special cases include the inviscid limit of viscous models and the a — > 
limit in a models. Next we show existence of a global attractor for the general model, 
and then give estimates for the dimension of the global attractor and the number of de- 
grees of freedom in terms of a generalized Grashof number. We then establish some 
results on determining operators for the two distinct subfamilies of dissipative and non- 
dissipative models. We finish by deriving some new length-scale estimates in terms of 
the Reynolds number, which allows for recasting the Grashof number-based results into 
analogous statements involving the Reynolds number. In addition to recovering most 
of the existing results on existence, regularity, uniqueness, stability, attractor existence 
and dimension, and determining operators for the well-known specific members of this 
family of regularized Navier-Stokes and MHD models, the framework we develop also 
makes possible a number of new results for all models in the general family, including 
some new results for several of the well-studied models. Analyzing the more abstract 
generalized model allows for a simpler analysis that helps bring out the core common 
structure of the various regularized Navier-Stokes and magnetohydrodynamics models, 
and also helps clarify the common features of many of the existing and new results. 
To make the paper reasonably self-contained, we include supporting material on spaces 
involving time, Sobolev spaces, and Gronwall-type inequalities. 
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1. Introduction 

The mathematical theory for global existence and regularity of solutions to the three- 
dimensional Navier-Stokes equations (3D NSE) is considered one of the most challeng- 
ing unsolved mathematical problems of our time [16]. It is also well-known that direct 
numerical simulation of NSE for many physical applications with high Reynolds number 
flows is intractable even using state-of-the-art numerical methods on the most advanced 
supercomputers available. Over the last three decades, researchers have developed tur- 
bulence models as an attempt to side-step this simulation barrier; the aim of turbulence 
models is to capture the large, energetic eddies without having to compute the smallest 
dynamically relevant eddies, by instead modeling the effects of small eddies in terms of 
the large scales in both NSE and magnetohydrodynamics (MHD) flows. 

In 1998, the globally well-posed 3D Navier-Stokes-a (NS-a) equations (also known as 
the viscous Camassa-Holm equations and Lagrangian averaged Navier-Stokes-a model) 
was proposed as a sub-grid scale turbulence model [8, 9, 10, 19, 20, 34, 51, 56] (see 
also [64] for n-dimensional viscous Camassa-Holm equations in the whole space). The 
inviscid and unforced version of 3D NS-a was introduced in [34] based on Hamilton's 
variational principle subject to the incompressibility constraint div u = (see also [52]). 
By adding the correct viscous term and the forcing / in an ad hoc fashion, the authors 
in [8, 9, 10] and [20] obtain the NS-a system. In references [8, 9, 10] it was found that 
the analytical steady state solutions for the 3D NS-a model compared well with averaged 
experimental data from turbulent flows in channels and pipes for a wide range of large 
Reynolds numbers. It was this fact which led the authors of [8, 9, 10] to suggest that 
the NS-a model be used as a closure model for the Reynolds-averaged equations. Since 
then, it has been found that there is in fact a whole family of globally well-posed 'a'- 
models which yield similarly successful comparisons with empirical data; among these 
are the Clark-a model [4], the Leray-a model [12], the modified Leray-a model [38], 
and the simplified Bardina model [5, 46]. 
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In addition to the early success of the a-models mentioned above, the validity of the 
original a-model (namely, the NS-a model) as a sub-grid scale turbulence model was 
also tested numerically in [11, 56]. In the numerical simulation of the 3D NS-a model, 
the authors of [11, 27, 28, 56] showed that the large scales of motion bigger than a (the 
length scale associated with the width of the filter which regularizes the velocity field) 
in a turbulent flow are captured (see also [49, 50] for the 2D case and [7] for the rate of 
convergence of the 2D a-models to NSE). For scales of motion smaller than the length 
scale a, the energy spectra decays faster in comparison to that of NSE. This numerical 
observation has been justified analytically in [19]. In direct numerical simulation, the fast 
decay of the energy spectra for scales of motion smaller than the supplied filter length 
represents reduced grid requirements in simulating a flow. The numerical study of [11] 
gives the same results. The same results hold as well in the study of the Leray-a model 
in [12, 27]. 

The NS-a turbulence model has also been implemented in a primitive equation ocean 
model (see [32, 33, 60]). Their simulations with the NS-a in an idealized channel domain 
was shown to produce statistics which resembles doubling of resolution. For other ap- 
plications of a regularization techniques, see [1] for application to the quasi-geostrophic 
equations, [45] for application to Birkhoff-Rott approximation dynamics of vortex sheets 
of the 2D Euler equations, and [48, 54, 55] for applications to incompressible magneto- 
hydrodynamics equations. In [6], an a-regularized nonlinear Schrodinger equation was 
proposed for the purpose of a numerical regularization that is hoped to shed some light 
on the profile of the blow-up solutions to the nonlinear Schrodinger equations. Also, 
in [2], the authors extend the derivation of the inviscid version of NS-a (called Euler-a, 
also known as Lagrangian averaged Euler-a) to barotropic compressible flows. 

Perhaps the newest addition to the family of a turbulence model is the Navier-Stokes- 
Voight (NSV) equations proposed in [5], which turn out to also model the dynamics of 
Kelvin- Voight viscoelastic incompressible fluids as introduced in [58, 59]. The statistical 
properties of 3D NSV have been studied in [47]. The long-term asymptotic behavior of 
solutions is studied in [42, 43]. In [18], the NSV was used in the context of image in- 
painting. The numerical study of NSV in [18] suggests that the NSV, in comparison with 
NSE, can provide a more efficient numerical process when automating the inpainting 
procedure for certain classes of images. It is worthwhile to note that the inviscid NSV 
coincide with the inviscid simplified Bardina model which is shown in [5] to be globally 
well-posed. This new regularization technique for Euler equations prevents the risk of 
damping too much energy in the small scales which could lead to unrealistic numerical 
results. 

As a representative of the more general model considered in this paper (described in 
detail in Section 2), we consider first the following somewhat simpler constrained initial 
value problem on an 3-dimensional flat torus T 3 : 

d t u + Au + (Mu ■ V)(Nu) + xV(Mm) t • (Nu) + Vp = f(x), 

V-u = 0, (1.1) 

u(0) = M , 

where A, M, and N are bounded linear operators having certain mapping properties, and 
where \ is either 1 or 0. As in prior work on regularized models of the Navier-Stokes, 
and Euler equations, we employ a single real parameter 9 to control the strength of the 
dissipation operator A. We then introduce two parameters which control the degree of 
smoothing in the operators M and N, namely 9\ and 9 2 , respectively, when x = 0, and 
6*2 and 9\, respectively, when x = 1. Some examples of operators A, M, and N which 
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satisfy the mapping assumptions we will need in this paper are 

A = (-A) e , M = (I -a 2 A)- e \ N = (I -a 2 A)- 62 , (1.2) 

for fixed positive real number a and for specific choices of the real parameters 8, 6\, and 
$2. However, we emphasize that the abstract mapping assumptions we employ are more 
general, and as a result do not require any specific form of the parameterizations of A, 
M, and N; this abstraction allows (1.1) to recover most of the existing regularization 
models that have been previously studied, as well as to represent a much larger three- 
parameter family of models that have not been explicitly studied in detail. As a result, the 
system in (1.1) includes the Navier-Stokes equations and the various previously studied 
a turbulence models as special cases, namely, the Navier-Stokes - a model [8, 9, 10, 19, 
20, 34, 56], Leray-a model [12], modified Leray-a model (ML -a) [38], simplified 
Bardina model (SBM) [5], Navier-Stokes- Voight (NSV) model [42, 43, 58, 59], and the 
Navier-Stokes -a- like (NS - a- like) models [57]. For clarity, some of the specific well- 
known regularization models recovered by (1.1) for particular choices of the operators 
A, M, N and x are listed in Table 1 . 



TABLE 1. Some special cases of the model (1.1) with a > 0, and with S — (I — 

aA)- 1 and Sg 2 = [I + (-aA) 8 *}- 1 . 
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Our main goal in this paper is to develop well-posedness and long-time dynamics 
results for the entire three-parameter family of models, and then subsequently recover 
the existing results of this type for the specific regularization models that have been 
previously studied. Along these lines, we first establish a number of results for the entire 
three -parameter family, including results on existence, regularity, uniqueness, stability, 
linear and nonlinear perturbations (with the inviscid and a — > limits as special cases), 
existence and finite dimensionality of global attractors, and bounds on the number of 
determining degrees of freedom. Elaborating on the latter a bit more, for 9 > 0, we 
derived a lower bound for the number of degrees of freedom m given by m > G n ^ e , 
where G is the Grashof number and n is the spatial dimension. A lower bound for the 
nondissipative case is also established. These results give necessary and/or sufficient 
conditions on the ranges of the three parameters for dissipation and smoothing in order 
to obtain each result, and we indicate where appropriate which particular regularization 
models are covered in the allowable parameter ranges for each result. In the final section 
of the paper, we develop some tools for relating the Grashof number-based results to 
analogous statements involving the Reynolds number. Analyzing a generalized model 
based on abstract mapping properties of the principal operators A, M, and iV allows 
for a simpler analysis that helps bring out the core common structure of the various 
regularized NSE (as well as regularized magnetohydrodynamics) models, and also helps 
clarify the common features of many of the existing and new results. 

In [57], a two-parameter family of models was studied, corresponding to a subset of 
those studied here, which we will call here the NS-a-like models. In order to describe this 
subset of models, let 9 and 6*2 be two nonnegative parameters, and consider the following 
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system on T 3 : 



d t u + (—A) it + (Mu ■ V)u + V(Mu) T -u + Vp 

V ■ u 
V ■ (Mu) 
«(0) 











(1.3) 



where M = (I + (— c^A) 6 * 2 ) -1 . This family of NS-a-like model equations interpolates 
between incompressible hyper-dissipative equations and the NS-a models when varying 
the two nonnegative parameters 8 and 8 2 . This is a special case of (1.1) with 81 = 
and x — h with the degree of dissipation controlled by the parameter 8 and the degree 
of nonlinearity controlled by only one parameter 8 2 . In this particular case, the NSE are 
obtained when 8 = 1 and 8 2 = 0, while the NS-a model is obtained when 8 = 8 2 = 1. 
In [57], sufficient conditions on the relationship between 8 and 8 2 are established to 
guarantee global well-posedness and global regularity of solutions. Our result here can 
be viewed as generalizing the global well-posedness and regularity results in [57] to 
a larger three -parameter family using a more abstract framework that does not impose a 
specific form for the parameterizations, and then also establishing a number of additional 
new results for the larger three-parameter family, including results on stability, linear and 
nonlinear perturbations, existence and dimension of global attractors, and on determining 
operator bounds. 

As a subset of the results mentioned above, we list some of the new results that we have 
obtained for the family of a models as a special case of the more generalized equation 
(1.1). As far as we know, these results have not been previously established in the liter- 
ature. The global existence and uniqueness of solutions for the inviscid a sub-grid scale 
turbulence models has been established only for the SBM [5]. Here, as a consequence 
of a more general result, we have established the global existence of a weak solution to 
the inviscid Leray-a-model of turbulence. In [20], the convergence of weak solutions of 
the NS-a to a weak solution of the NSE as the parameter a— >0 was established. Here we 
have established this convergence result as well for the NS-a-like equations. In addition, 
we have established for the NS-a-like equations, the existence and finite dimensionality 
of its global attractor, and determining operator bounds. In the case of Leray-a, ML-a, 
and SBM, the results on determining operator bounds also appear to be new. 

It is important to note that the general framework here allows for the development of 
results for certain (regularized or un-regularized) magnetohydrodynamics (MHD) mod- 
els. The basic MHD system has the form 



where the unknowns are the velocity field u, the magnetic field h, and the pressure it, 
and where v > and i] > denote the constant kinematic viscosity and constant mag- 
netic diffusivity, respectively. Our global well-posedness results include for example a 
particular regularized MHD model 



d t u — uAu + (u ■ V)u 
d t h - r]Ah + (u ■ V)h 
V 



- (h-V)h 

- (h ■ V)u 
h = V -u 



Vvr 




0. 

0. 



(1.4) 



d t u - uAu + (Mu ■ V)« - (Mh ■ V)h = Vvr - ^V|/i| 2 , 
d t h - r)Ah + (Mu ■ V)h - (Mh ■ V)u = 0, 



(1.5) 
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supplemented with several divergence-free and boundary conditions, where M = (I — 
a 2 A)' 1 . Note that the system (1.5) is different from the 3D Leray-a-MHD model pro- 
posed in [48], where global well-posedness result is still open as in the case of the original 
MHD equations. It is also different in nature from the MHD model proposed in [63]. For 
the 3D Leray-a-MHD model proposed in [48] the magnetic field h is not regularized. 
Another regularized model whose global well-posedness result is covered here is the 
following modified version of the MHD-a model proposed in [48], 



supplemented with several divergence-free and boundary conditions. Again, the above 
system is different from the original version of the MHD-q system proposed in [48], 
which does not have a regularization on the magnetic field h. For the original MHD-a 
system, global well-posedness is established in [48]. Here we would like to stress that 
our current framework is best suited to MHD models where the velocity field u and the 
magnetic field h are treated on an equal footing as far as the regularization is concerned, 
so we cannot obtain sharpest results in our framework for the systems like the ones 
proposed in [48]. However, our framework requires only minor modification to include 
these models; the function spaces become product spaces, and the principal dissipation 
and smoothing operators become block operators on these product spaces, typically with 
block diagonal form. It is worthwhile to note here that filtering the magnetic field, as is 
done in [54, 55], can be interpreted as introducing hyperviscosity for the filtered magnetic 
field. This observation was first introduced in [48], and smoothing the magnetic field was 
thought to yield unnecessary extra dissipation added to the system. 

Since the a models of turbulence were intended as a basis for regularizing numerical 
schemes for simulating turbulence, it is important to verify whether the ad hoc smoothed 
equations inherit some of the original properties of the 3D MHD equations. In particular, 
one would like to see if the ideal (v — r\ — 0) quadratic invariants of the smoothed MHD 
system can be identified with the ideal invariants of the original 3D MHD system under 
suitable boundary conditions. There are three ideal quadratic invariants in 3D MHD 
(e.g. under rectangular periodic boundary conditions or in the whole space), namely, the 
energy E = | J n \u(x)\ 2 + \h(x)\ 2 dx, the cross helicity h c = | J n u(x) ■ h(x) dx, 
and the magnetic helicity h M = | J Q a(x) ■ h(x) dx, where a(x) is the vector potential 
so that h(x) = V x a(x). In the case of the 3D MHD-a equations in [48], the three 
corresponding ideal invariants are the energy E a = \ J n u(x) ■ Mu(x) + |/i(a;)| 2 dx 
(which reduce, as a — » 0, to the conserved energy of the 3D MHD equations), the 
cross helicity — | J Q u(x) ■ h(x) dx, and the magnetic helicity = \ J n a(x) ■ 
h(x) dx. For our system in (1.5) the corresponding ideal invariants are the energy E a = 
| f n \u(x)\ 2 + \h(x)\ 2 dx, and the cross helicity = \ J n u(x) ■ h(x) dx. Currently, we 
are unable to identify a conserved quantity in the ideal version of (1 .5) which correspond 
to the magnetic helicity in the 3D MHD system. Similar problems arise in the 3D MHD- 
Leray-a equations introduced in [48]. When both the magnetic field and the velocity 
field are filtered as it is done in [54], the corresponding ideal quadratic invariants are 
E a = | J n u(x)-Mu(x)+h(x)-Mh(x) dx, the cross helicity h% = | J n u(x)-Mh(x) dx, 
and the magnetic helicity = \ j n Ma(x) ■ Mh(x) dx. 

The remainder of the paper is structured as follows. In §2, we establish our notation 
and give some basic preliminary results for the operators appearing in the general reg- 
ularized model. In §3, we build some well-posedness results for the general model; in 
particular, we establish existence results (§3.1), regularity results (§3.3), and uniqueness 



\u - uAu + (Mu ■ V)tt + V(Mu) T ■ u - (Mh ■ V)h = Vvr 
d t h - rjAh + (Mu ■ V)A - (Mh ■ V)u = 0, 




(1.6) 
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and stability results (§3.2). In §4 we establish some results for singular perturbations, 
which as special cases include the inviscid limit of viscous models and the a — ► 
limit in a models; this involves a separate analysis of the linear (§4.1) and nonlinear 
(§4.2) terms. These well-posedness and perturbation results are based on energy meth- 
ods. In §5, we show existence of a global attractor for the general model by dissipation 
arguments (§5.1), and then by employing the classical approach from [62], give esti- 
mates for the dimension of the global attractor (§5.2). In §6, we establish some results 
on determining operators for the two distinct subfamilies of dissipative (§6.1) and non- 
dissipative (§6.2) models, with the help of certain generalizations of the techniques used 
e.g., in [21, 41, 37, 43]. Since the results in §6 are (naturally) given in terms of the 
generalized Grashof number, we finish in §7 by developing some new results on length- 
scale estimates in terms of the Reynolds number, which allows for relating the Grashof 
number-based results in the paper to analogous statements involving the Reynolds num- 
ber. 

To make the paper reasonably self-contained, in Appendix A we develop some sup- 
porting material on Gronwall-type inequalities (Appendix A.l), spaces involving time 
(Appendix A.2), and Sobolev spaces (Appendix A. 3). 

2. Notation and preliminary material 

Let be an n-dimensional smooth compact manifold with or without boundary and 
equipped with a volume form, and let E — > f2 be a vector bundle over £1 equipped 
with a Riemannian metric. With C°°(E) denoting the space of smooth sections of E, let 
V C C°° (E) be a linear subspace, let A : V— ► V be a linear operator, and let 73 : V x V^V 
be a bilinear map. At this point V is conceived to be an arbitrary linear subspace of 
C°°(E); however, later on we will impose restrictions on V implicitly through various 
conditions on certain operators such as A. Assuming that the initial data u £ V, and 
forcing term / £ C°°(0, T; V) with T > 0, consider the following equation 

d t u + Au + B(u, u) = f, 

U[0) = M , 

on the time interval [0, T]. Bearing in mind the model (1.1), we are mainly interested in 
bilinear maps of the form 

B(v,w) = B(Mv,Nw), (2.2) 

where M and N are linear operators in V that are in some sense regularizing and are 
relatively flexible, and B is a bilinear map fixing the underlying nonlinear structure of 
the equation. In the following, let P : C°°(E)^V be the L 2 -orthogonal projector onto 
V. 

Example 2.1. a) When Q is a closed Riemannian manifold, and E = TQ the tangent 
bundle, an example ofV is V pcr C {u £ C°°(TQ) : divu = 0}, a subspace of the 
divergence-free functions. The space of periodic functions with vanishing mean on a 
torus T n is a special case of this example. In this case, one typically has A = (—A) 9 , 
M = (I - a 2 A)- e \ and N = {I - a 2 A)^ 2 . 

b) When Q is a compact Riemannian manifold with boundary, and again E = TQ the 
tangent bundle, a typical example ofV is Vh om = { u £ C^(TQ) : divu = 0} the space 
of compactly supported divergence-free functions. In this case, we keep the operators 
A = (—PA) 8 , M = (I -a 2 PA)- 01 , andN = (I -a 2 PA)- 02 , in mind as the operators 
that one would typically consider. 
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c) In either of the above two examples, one can consider as V the product spaces 
Vper x V pcr and Vh om x Vhom> which are encountered, e.g., in magnetohydrodynamics 
models, cf. Example 2.2 below. For the operators A, M, and N, we would have the 
corresponding block operators on the product space V for the above two examples, acting 
diagonally. 

Example 2.2. a) In a) or b) of the Example 2.1 above, the bilinear map B can be taken 
to be 

B x {v,w) = P[(v ■ V)w], (2.3) 

which corresponds to the models with x = as discussed in §i. 

b) Again in a) or b) of Example 2.1 above, B can be taken to be 

B 2 (v, w) = P[(w ■ V)« + (V/)u], (2.4) 

which corresponds to the models with x — 1 as discussed in §i. 

c) An example of B that cannot be written in the form (2.2) is the bilinear map for the 
Clark-a model, which is 

B c (v,w) = B 1 (Nv,w) + B 1 (w,Nv) - B^Nv, Nw) - aB 1 (\/ j Nv,VjNw), 

where N = (I — a^PA)^ 1 , and where the Einstein summation convention is assumed. 
Note that for this bilinear map, one only has (B c (v,v), Nv) = for any v e V per or 
v G Vhom. in contrast to the examples B\ and B 2 which are well-known to have stronger 
antisymmetry properties. 

d) The MHD system (1.4) has the bilinear map 

B 3 (v,w) = (B^vt, wx) - Bi(v 2 ,w 2 ), Bi(v u w 2 ) - Bx^wx)) , (2.5) 

where v = (v i, v 2 ) and w = (wx,W2) are elements of either V per x V per or Vhom x Vh om - 
The bilinear map for (1.5) is a special case of B 3 (v, w) = B 3 (Mv, Nw), with M and 
N having the form M = diag(Mi, Mi) and N = I. Another special case of B 3 is the 
bilinear form for the Leray-a-MHD model proposed in [48], where M and N have the 
form M = diag(Mi, I) and N = I. 

e) The MHD-a model (1.6) has the bilinear map of the form 

B 4 (v,w)= (B 2 (v 1 ,w 1 ) - Bi(v2,W2),B 1 (v 1 ,w 2 ) - B\ (v 2 , Wi)) • (2.6) 

The bilinear map for (1.6) is a special case of B^v, w) = B±(Mv, Nw), with M and 
N having the form M = diag(M l5 Mi) and N = I. Another special case of B^ is the 
bilinear form for the MHD-a model proposed in [48], where M and N have the form 
M = diag(Mx, I) and N = I. 

f) More generally, one can consider the bilinear maps of the form 

B B (v,w) = B ijjyk (v,w) = (5i(t»i,wi) - Bj{v 2 ,W2),Bk{vi,w 2 ) - Bj(v 2 , wi)) , 

where i,j, k e {1, 2}. This class includes the above examples d) and e). 

To refer to the above examples later on, let us introduce the shorthand notation: 

Bi(v,w) = Bi(Mv,Nw), i = l,...,5. (2.7) 

As usual, we will study equation (2.1) by extending it to function spaces that have 
weaker differentiability properties. To this end, we interpret the equation (2.1) in distri- 
bution sense, and need to continuously extend A and B to appropriate smoothness spaces. 
Namely, we employ the spaces V s = closes V, which will informally be called Sobolev 
spaces in the following. The pair of spaces V s and V~ s are equipped with the duality 
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pairing (-, ■), that is, the continuous extension of the L 2 -inner product on V°. More- 
over, we assume that there is a self-adjoint positive operator A such that A s : V s — > V° 
is an isometry for any s E E. For arbitrary real s, assume that A, M, and A" can be 
continuously extended so that 

A:V S ^V S ~ 29 , M :V s ^V s+2dl , and N:V s ^V s+2e \ (2.8) 

are bounded operators. (Again, the assumptions we will need for A, M, and A" are more 
general, and do not require this particular form of the parameterization.) We will assume 
9 > 0; however, there will not be any a priori sign restriction on 9\ and 9 2 . We remark 
that s in (2.8) is assumed to be arbitrary for the purposes of the discussion in this section, 
and that it is of course sufficient to assume (2.8) for a limited range of s for most of the 
results in this paper. 

Remark 2.3. Note that in this framework the best value for 9\ is 9\ = for both the 
Leray-a-MHD model and the MHD-a model as proposed in [48], since those models 
have M = (Mi, I), cf. Example 2.2 d) and e). It is possible to refine our analysis by 
considering spaces such as V s x V r instead of V s x V s . 

We assume that A and A" are both self-adjoint, and coercive in the sense that for (3 E E, 

(Av,A 2/3 v) > c A \\v\\ 2 e+ p - C A \\v\\% v E V 9 ^, (2.9) 

with c A = c A {j3) > 0, and C A = C A ((3) > 0, and that 

(Nv,v) > c N \\v\\ 2 _ 92 , vEV~ e \ (2.10) 

with cat > 0. We also assume that (2.9) is valid for (3 = —9 2 with A 2/3 replaced by 
N. Note that if 9 = 0, (2.9) is strictly speaking not coercivity and follows from the 
boundedness of A, and note also that (2.10) implies the invertibility of N. For clarity, we 
list in Table 2 the corresponding values of the parameters and bilinear maps discussed 
above for special cases listed in Table 1 . 

TABLE 2 . Values of the parameters 8, 9\ and 8 2 , and the particular form of the bilinear 
map B for some special cases of the model (2.1). (The bilinear maps B\ and B 2 are as 
in (2.3)-(2.4), (2.7).) Each of the NSE models has a corresponding MHD analogue. 



Model 


NSE 


Leray-a 


ML-o 


SBM 


NSV 


NS-a 


NS-a-like 


6 


1 


1 


1 


1 





1 


9 







1 





1 


1 
















1 


1 


1 


1 


2 


B 


#1 


Bi 


Bi 


B x 


Bi 


B 2 


B 2 



We denote the trilinear form b(u, v, w) = (B(u, v),w), and similarly the forms b, bi, 
and bi, i = 1, . . . , 5, following Example 2.2 and (2.7). We consider the following weak 
formulation of the equation (2.1): Find u E Ll oc (0, T; V s ) for some s such that 

/ (-(u(t),w(t)) + (Au{t),w{t)) + b{u{t),u{t),w{t)) - (f{t),w{t)) )dt = 
Jo 

U(0) = M , 

(2.11) 

for any w E C^(0, T; V). Here the dot over a variable denotes the time derivative. 

The left hand side of the first equation in (2.11) is well defined if u E L 2 (0, T; V s ), 
f E L l {f),T]V s '), and b : V s x V s x V~<-*R is bounded for some s,s',7 E E. The 
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second equation makes sense if u G C(0, e; V s ) for some e > and s6R. However, 
the following lemma shows that the latter condition is implied by the first equation. 

Lemma 2.4. Let X C L 2 (0, T; V s ) be the set of functions that satisfy the first equation 
in (2.11). Let f G L\0, T; V s - 29 ), and let b : V s x V s x V 2e - S ^R be bounded. Then 
we have X C C(0, T; V s ~ 2e ). 

Proof. IfueX, we have u G -^((^T; V s ~ 2d ); and so [61, Lemma 3.1.1] implies that 

u G C(0,T;V s - 29 ). □ 

In concluding the preliminary material in this section, we state the following result on 
the trilinear forms bi, i = 1, . . . , 5, which are the main concrete examples for the trilinear 
form b. Recall that b(u, v, w) = b(Mu, Nv, w). 

Proposition 2.5. a) The trilinear forms \, % = 1, . . . , 5, are antisymmetric in their sec- 
ond and third variables: 

b~i(w,v,v) = 0, w,vEV, 2 = 1,..., 5, (2.12) 

where V is one of the appropriate spaces in Example 2.1. 

b) The trilinear form bi : V ai x V 72 x V 73 — >R is bounded provided that 

a 1 + a 2 + a 3 >^, (2.13) 

and that for some k G {0, 1}, 

0"2+ (T 3 > 1, CTi + (T3>A;, (7i + (7 2 > 1 - k. (2.14) 

//■^e above three conditions are satisfied, and if (Tj is a non-positive integer for some 
i G {1, 2, 3}, ?/ze inequality in (2.13) can be replaced by the non-strict version of 
the inequality. The non-strict inequality is also allowed if for some k G {0, 1}, 

0"! >0, a 2 > k, a 3 >l-k. (2.15) 

c) The trilinear form l>2 '■ V ai x V a ' 2 x V a3 — >M is bounded provided that (2.13) holds 
and in addition that 

&2 + ^3 > 0, cri+<7 3 >l, a 1 + a 2 >l. (2.16) 

If the above three conditions are satisfied, and if <j{ is a non-positive integer for some 
i G {1, 2, 3}, then the inequality in (2.13) can be replaced by the non-strict version of 
the inequality. The non-strict inequality is also allowed if 

o-i >1, 02 >0, a 3 >0. (2.17) 

d) The trilinear form 63 : V ai x V 72 x V 73 — >M is bounded under the same conditions 
on <7i, (J2, and 03 that are given in b). 

e) The trilinear forms 64, 65 : V U1 x V u ' 2 x V^ 73 — *>IR are bounded provided that (2.13) 
holds and in addition that 

02 + 03 >L 0i + 3 >l, 0i + 02>L (2.18) 

If the above three conditions are satisfied, and if Oi is a non-positive integer for some 
i G {1, 2, 3}, then the inequality in (2.13) can be replaced by the non-strict version of 
the inequality. The non-strict inequality is also allowed if for some k G {0, 1}, 

01 > 1, 02 > k, 0-3 >1- k. (2.19) 
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Proof. The antisymmetricity of bi, and b 2 is well known, and the boundedness of Bi is 
immediate from Lemma A.3. The antisymmetricity of 6 5 (which a fortiori implies that 
of 63 and 64) can be seen from 

b 5 (w, v,v) = b i (w 1 ,v 1 ,v 1 ) - b j (w 2 ,v 2 ,v 1 ) + b k {wi,v 2 ,v 2 ) - b j (w 2 ,v 1 ,v 2 ), 

where i, j, k G {1, 2}. Part b) is proven by applying Lemma A. 3 to each of the following 
two representations 



Similar to the Leray theory of NSE, we begin the development of a solution theory for 
the general 3-parameter family of regularized Navier-Stokes and MHD models with clear 
energy estimates that will be used to establish existence and regularity results, and un- 
der appropriate assumptions show uniqueness and stability. To reinforce which existing 
results we recover in this general unified analysis, we give the corresponding simplified 
results that have been established previously in the literature for the special cases listed 
in Table 1, at the end of the proof of every theorem. 

3.1. Existence. In this subsection, we establish sufficient conditions for the existence of 
weak solutions to the problem (2.1 1). By a weak solution, we mean a solution satisfying 

u G L 2 (0, T; V 9 ~ 9 ' 2 ) and u G -^ 1 (0, T; V 1 ) for some 76M. 

Theorem 3.1. a) Let the following conditions hold. 

i) b : V ai x V 72 x V-tR is bounded for some <t { G [-9 2 , 9 - 9 2 \, i = 1,2, and 



7 g [9 + 9 2 , 00) n (#2,00); 

ii) b(v,v,Nv) = for any v G V ' 02 ; 

iii) b : V ai x V a ' 2 x V^—tR is bounded for some Oi < 9 — 9 2 , % = 1, 2, and 7 > 7; 

iv) u G V- 02 , and f G L 2 (0, T; V' - 02 ), T > 0. 

Then, there exists a solution u G L°°(0, T; V~ 02 ) n L 2 (0, T; V ' 02 ) to (2.11) satisfying 



b) With some (3 > —9 2 , let the following conditions hold. 

i) b : V 13 x V 13 x V e ~ p ^R is bounded; 

ii) M e V?, and f G L 2 (0, T; V~ 0+l3 ), T > 0. 

iii) b : V a x V 7 x V 7 ^IR is bounded for some o < 9 + (3, and 7 > 9 — j3; 

Then, there exist T*(u ) = T*(\\u \\p) > and a local solution u G L°°(0,T*; V&) n 
L 2 (0, T*; V r6l+/3 ) to equation (2.11). 

Remark 3.2. a) Let 9 + #1 > |. Then from Proposition 2.5 the trilinear forms b\ and 63 
/m/^Z/ ?/ze hypotheses of a) for < 9 -9 2 -l with -7 < min{26> + 26»i - 6» - 02 + 
26 1 !, # + ^2 — 1}. A^ofe m particular this gives the global existence of a weak solution 
for the inviscid Leray-a model. As far as we know this result has not been reported 
previously. 



bx{u,v,w) = (u l ViV k ,w k ) = (u i v k ,V i w k ) . 
Part c) is proven by applying Lemma A.3 to 



b 2 {u, v, w) = <^ViM fc + u l V k v h w k ) = (v i V i u k + V^M^i) - ^V fe W;, w k ) 
= (v i V i u k -ViV k Ui,w k ) . 
To complete the proof, d) and e) follow from parts b) and c). 



□ 
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b) Let 9 + 20x > k, 9 + 20 2 > 1, (3 > ^ - 9 - 2{9 X + 2 ), and (3 > 1=* - 9 X - 9 2 , 
for some k G {0, 1}. Then the trilinear forms b x and 63 satisfy the hypothesis ofb). 

Remark 3.3. a) Let 9 + 9\ > |. Then the trilinear form b 2 fulfills the hypotheses of a) 
for -7 < - 2 - 1 with -7 < min{20 + 20i - =±*, - 2 + 20i - 1, + 2 }. 

b) The trilinear form b 2 satisfies the hypothesis of b) for (3 > lt ^ L — 9 — 2(9 1 + 2 ) with 
/3 > I - 01 - 9 2 provided 9 + 29 x >land9 + 20 2 > 0. 

Remark 3.4. a) Let 9 + 9\ > \. Then the trilinear forms 64 and b 5 fulfill the hypotheses 
of a) for -7 < - 2 - 1 with -7 < min{20 + 20! - *±*, - 2 + 20 x - 1, + 2 - 1}. 

b) Let 9 + 20 x > 1, + 20 2 > 1, (3 > - - 2(0i + 2 ), and (3 > ± - 9 X - 2 . 
TTzen f/ze trilinear forms b\ and 63 satisfy the hypothesis ofb). 

Proof of Theorem 3.1. Let {Vm : m G N} C V 6 ^ 6 * 2 be a sequence of finite dimensional 
subspaces of V 9 ~ 92 such that 

(1) V m C V m+1 for all m G N; 

(2) Li m( z N V m is dense in V e ~ d2 ; 

(3) For m G N, with W TO = NV m C ^ e+e2 , the projector P m : K e " e2 ^\/ m defined 
by 

(P m v, w m ) = (v, w rn ) , w m G W m , i> G V, 
is uniformly bounded as a map V^ 7 ^V~ 7 . 
Such a sequence can be constructed e.g., by using the eigenfunctions of the isometry 

Consider the problem of finding u m G C 1 (0, T; V m ) such that for all w m G W m 

(u m , w m ) + (Au m , w m ) + b(u ) = (f,W m ), 

(u m (0),w m ) = (u ,w m ). 

Upon choosing a basis for V m , the above becomes an initial value problem for a system 
of ODE's, and moreover since iV is invertible by (2.10), the standard ODE theory gives 
a unique local-in-time solution. Furthermore, this solution is global if its norm is finite 
at any finite time instance. 

The second equality in (3.1) gives 

CN\\u m m\-82 < (u m (0),Nu m (0)) = (u(0),Nu m (0)) < ||u(0)||_ 9 JiVu m (0)|| e2 , 
so that 

R»(0)||-fl.< im ' d2 '' e2 \\u(0)\\-e 2 . (3-2) 

c N 

Now in the first equality of (3.1), taking w m = Nu m , and using the condition ii) on b, we 
get 

^(u m , Nu m ) + 2(Au m , Nu m ) = 2(f m , Nu m ) 

^^ll/ll-o-^ +£11^11^11^11^, 
for any e > 0. Since by choosing e > small enough we can ensure 

2( i AlL m , Nu m ) + s\\ iV||_0 2 .0 2 ll^mllfl-^ ~ ll^mll— 2 > 

by Gronwall's inequality we have 

\\u m (t)\\ 2 -02 £ (lK(0)||% 2 + f ll/ll V J e c \ (3.4) 
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for some C E R. For any fixed T > 0, this gives u m E L°°(0, T; V~ 02 ) with uniformly 
(in m) bounded norm. Moreover, integrating (3.3), and taking into account (3.4), we 
infer 

/ (Au m ,Nu m ) <ip(t), *e[0,oo), (3.5) 
Jo 

where ij) : [0, oo)— >(0, oo) is a continuous function. If 9 > 0, by the coerciveness of 
A, the above bound implies u m E L 2 (0, T; V ~ 02 ) with uniformly bounded norm. So in 
any case, u m is uniformly bounded in L°°(0, T; V~ 02 ) R L 2 (0, T; V ~ 02 ), and passing to 
a subsequence, there exists u E L°°(0, T; V~ 02 ) H L 2 (0, T; V^~ 92 ) such that 

u m -> u weak-star in L°°(0, T; V" 02 ), 

u m — > ti weakly in L 2 (0,T; y 9-62 ). D ' 6) 

For passing to the limit m— >oo in (3.1) we shall need a strong convergence result, 
which is obtained by a compactness argument. We proceed by deriving a bound on 
u m = ^. Note that (3.1) can be written as 

Um PmAu m ~\- P m B(u m ,U n i) P m f y 

Um(0) = P m u(0). 

Therefore 

||«m||-7 < C\\u m \\ e _ e2 + \\P m B(u m ,u m )\\^ + ||P m /||_ 7 =: I x + I 2 + 7 3 . (3.8) 
By the boundedness of B, we have 

■^2 W^m || cri ||w m ||o- 2 . (3.9) 

If 6* = 0, then the norms in the right hand side are the V _6,2 -norm which is uniformly 
bounded. If 9 > 0, since 



(3.7) 



^m||<Ti — ||w m ||_g 2 ||M m ||g_g 2 , Aj — n , i — 1)2, (3.10) 



|X-Ai tU, II A< \ _ a i + ^ 

9 

by the uniform boundedness of u m in L 00 ^^ 2 ) we have 

T < Wni || 2— Ai — A 2 || MA1+A2 < II, 11 A1+A2 /o |i\ 

J 2 ll u m||_e 2 \\ u "m\\0—0 2 r-u \\ a m\\o-e 2 " W- 1J J 

Hence, with A := Ai + A 2 = — ^ ^ 2 if 9 > 0, and with A = 1 if 9 = 0, we get 



Unroll jDP(y-7) ~ ll M w||^p(ye-8 2 ) + ll M m||^pA(ye-s 2 ) + ||/Hxp(y-fl-e2)- (3-12) 

The first term on the right-hand side is bounded uniformly when p < 2. The second term 
is bounded if p\ < 2, that is p < 2/ A. We conclude that w m is uniformly bounded in 

LP (0, T; V- - 02 ), with p = min{2, 2/A}. 

By employing Theorem A. 2, we can now improve (3.6) as follows. There exists u E 

C(0,T; V- - 02 ) nZ°°(0,T; V~ 02 ) nL 2 (0,T; V - 02 ) such that 

u m -> u weak- star in L°°(0,T; V -6 * 2 ), 

u m -f u weakly in L 2 (0, T; V 9 ^ 2 ), (3.13) 
u m u strongly in L 2 (0, T; V s ) for any s < 9 - 9 2 . 

Now we will show that this limit u indeed satisfies the equation (2.11). To this end, let 
w E C°°(0, T; V) be an arbitrary function with w(T) = 0, and let w rn E C\0, T; W m ) 
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be such that w m {T) = and w m ^w in C\0, T; V _J? ). We have 

(u m {t),w m {t))dt + / (Att m (t),w m (t))dt + / b{u m {t),u m {t),w m {t))dt 



= (u m (0),w m (0))+ (f(t),w m (t))dt. (3.14) 

</ o 

We would like to show that each term in the above equation converges to the correspond- 
ing term in 

(u(t),w(t))dt + [ (Au(t),w(t))dt+ [ b(u(t),u(t),w(t))dt 
Jo Jo 

= (u o ,w(0)) + [ (f(t),w(t))dt, (3.15) 
Jo 

which would imply that u satisfies (2. 11). Here we show this only for the nonlinear term. 
We have 

\b(u rn (t),u m (t),w m (t)) - b(u(t),u(t),w(t))\dt <I m + II m + III m , (3.16) 



where the terms I m , II m , and III m are defined below. Firstly, it holds that 



T 



< 



\b(u m (t),u m (t),w m (t) -w(t))\dt 
T (3.17) 

||Wm(*) ||*i ||Mm(*) |U 2 ||Wm(£) - w(t)\\jdt 



< ||Mm||L 2 (y s i)||w m ||L 2 (^2)||u' m — w\\ C (vy)- 
thus, we get linim^oo I m = since a fortiori o-i < 9 — 9 2 ,i = 1,2. For II m we have 

II m = / \b(u m (t) -u(t),u m (t),w(t))\dt 



< 



T 

\\u m (t) - u(t)\\ ai \\u m (t)\\ a2 \\w(t)\\ e+ e 2 dt 



(3.18) 



< \\u m — u\\ L 2^y<y 1 ' ) \\u m \\ L 2^ V a 2 ^\\w\\c(V"r), 

so linv^oo II m = since o\ < 9 — 9 2 and a 2 < 9 — 9 2 . Similarly, we have 

III m = [ \b{u(t),u m (t)-u{t),w{t))\dt 
Jo 

< ||«||l,2(yai)||u m - M|| i 2(y t T 2) ||u'||c(y7), 



(3.19) 



so linim^oo III m = since o\ < 9 — 9 2 and a 2 < 9 — 9 2 . 

For the proof of b), we choose the nested subspaces {V m : m E N} C V e+/3 and 
W m = ^ w V m C V e ~ p , in (3.1). Then substituting w m = A 2/3 w m (0) in the second 
equality in (3.1) gives 

Il«m(0)lli8 ^ hoh- 

Now in the first equality of (3.1), taking w m = A 2l3 u m , and using the boundedness of 6, 
we get 

^11 II 2 < II -fll 2 _, 1 1 II 2 , II ||4 
^ll^mll^ 11/ 11/3-0 I \\ U m\\f3 ' \\ u m\\p, 
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and thus 

for some T* > 0. The rest of the proof proceeds similarly to that of a). □ 

For clarity, the corresponding conditions and results of Theorem 3.1 above are listed 
in Table 3 below for the special standard model cases of the general three-parameter 
regularized model listed in Table 1. For the NS-a-like case in the table, the allowed 
values for (3 are (3 > § - 9 - 20 2 with (3 > \ - 9 2 provided 9 > 1 and 9 + 29 2 > 0. 

TABLE 3. Existence results for some special cases of the model (1.1). The table gives 
values of (a, b) and (j,p) for our recovered existence results for the standard models, 
giving existence of a solution u € L°°(V a ) n L 2 (V b ), with u <E L P (V"~ 7 ). (In the NSV 
case, e > 0, and in the NS-a-like case, 7 > max{-0 2 — 1/2, 6» 2 + 1 /2, n/2}.) The last 
row indicates the local existence of solutions u 6 L 00 (V r ' 3 ) n L 2 {V^ +e ). The result for 
each NSE model has a corresponding MHD analogue. 



Model 


NSE 


Leray-a 


ML-q 


SBM 


NSV 


NS-a 


NS-a-like 


Existence 
















a, b 


0, 1 


0, 1 


-1,0 


-1,0 


-1,-1 


-1,0 


— 02, — — 6 2 


7.P 


1 4 


1,2 


2,2 


2,2 


l + e,2 


2,2 


7,2 


Local 


P>\ 


> 


/3>-| 


P> -1 


/? > -1 


/?>-! 


/3 



3.2. Uniqueness and stability. Now we shall provide sufficient conditions for unique- 
ness and continuous dependence of on initial data for weak solutions of the general three- 
parameter family of regularized models. 

Theorem 3.5. Let /3 > -9 2 , and let u u u 2 G L°°(0,T; V p ) n L 2 {0,T; V p+e ) be two 
solutions of (2.11) with initial conditions Wi(0),tt2(0) G V 13 , respectively. 

a)Letb: V ai x V e ~ 92 x V n ->NL be bounded for some a x < 9 - 9 2 and a 2 < 9 + 9 2 
with <j\ + 02 < 9. Moreover, let b(v,w, Nw) = for any v G V ai and w G V U2 . Then 
we have 

|Mt) - u 2 {t)\\-e 2 < #)IM0) - « 2 (0)|U, t G [0,T], (3.20) 

where <p G C([0, T]) and 0(0) = 1. 

6) Le? 6 : V p x V 3 x V^-^R be bounded. Then we have 

\\ux{t) - u 2 (t)\\ p < 0(*)||«i(O) - w 2 (0)|| /3 , t G [0,T], (3.21) 

where <p G C([0, T]) and 0(0) = 1. 

Remark 3.6. The trilinear forms b\ and 63 satisfy the hypotheses of a) provided 9 + 61 > 
if* + 20 x > jfe, + 2 > §, 20 + 201 + 02 > ^ and 30 + 20! + 20 2 > 2-k,forsome 
k G {0, 1}. The forms bi and b 3 satisfy the hypothesis of b) for (3 > rk ^- — 2(9 1 + 2 ) — 
with (3 > i=* - (0! + 2 ) provided 20 2 + > 1 and 20 x + > fc, /or w fc G {0, 1} . 

Remark 3.7. The trilinear form b 2 satisfies the hypotheses of a) for 9 + 20i > l,9+9\ > 
|, + 2 > 0, 20 + 20x + 2 > s±2 an j 3^ + 20 x + 20 2 > 1. 77ze trilinear form b 2 
satisfies the hypothesis of b) for (3 > tk ^ — 2(0! + 2 ) — wM /3 > | — (0i + 2 ) provided 
20 2 + > and 20 x + > 1. 

Remark 3.8. Tne trilinear forms 64 an J 65 satisfy the hypotheses of a) provided 9 + X > 
i + 20i > 1, + 2 > \, 29 + 20x + 2 > =±?, ana" 30 + 20i + 20 2 > 2. The forms b 4 
and b 5 satisfy the hypothesis of b) for (3 > - 2(0 X + 2 ) - wd/z /3 > § - (0i + 02) 
provided 20 2 + > 1 and 20 x + > 1. 



16 



M. HOLST, E. LUNASIN, AND G. TSOGTGEREL 



Proof of Theorem 3.5. Let v — Ui — u 2 - Then subtracting the equations for U\ and u 2 we 
have 

(v, w) + (Av, w) + b(v, Mi, w) + b(u 2 ,v, w) = 0. (3.22) 
Taking w = Nv, we infer 
d„ „ 



dt l 



-8 2 



+ c\\v\L_n < C\\v 



e 2 ^ ^ imkilPilk 

< C1MI 2 ^ Al ~ A2|U ' l|Al+A2 

— II II— 02 



-02 11^11^2—2^2 
„.|| A1+A2 II „, || 

v\\o-e 2 \\ u i\\e-e 2 i 



where Ai 



0-1+6)2 



and A 5 



d 



By applying Young's inequality we get 



dt 11 ^ 
Now Gronwall's inequality gives 



< C||w||% 2 ||«i| 



-02- 



\Ht)\\U 2 < lb(O)||__ 02 exp / ckh^. 

./o 

The part b) is proven similarly, taking, e.g. w = (I — A)@v. 



(3.23) 

□ 



To clarify these results in the case of specific models, the corresponding conditions 
and results of Theorem 3.5 above are listed in the Table 4 below for the special case 
models listed in Table 1. 

TABLE 4. Uniqueness results for some special cases of the model (1.1). The table 
gives values of for our recovered uniqueness results for the standard models, where 
u e V fi and where u e L 00 (V^ /3 ) n L 2 (V l3+e ). (In the NS-a-like case, the requirement 
on (3 is that /3 > max{-6» 2 , 1/2 - 9 2 , 5/2 - 6 - 29 2 }). The result for each NSE model 
has a corresponding MHD analogue. 



Model 


NSE 


Leray-a 


ML-a 


SBM 


NSV 


NS-a 


NS-a-like 


Uniqueness 


P>\ 


p > 


(3>-\ 
or (3 = -1 


P>-1 


0>-l 


(3>-\ 
or 13= -1 


P 



3.3. Regularity. In this subsection, we develop a regularity result on weak solutions for 
the general family of regularized models. 

Theorem 3.9. Let u G L 2 (0, T; V e ~ 62 ) be a solution to (2.11), and with some (3 > —9 2 , 
let the following conditions hold. 

i) b : V a x V a x V ~ P -^R is bounded, where a = min{/3, 9 - 9 2 }; 
h) b(v, w, Nw) = Ofor any v, w G V; 
iii) w(0) G V p , andfe L 2 {0,T; V^ e ). 

Then we have 

u G L°°(0, T; V 13 ) fl L 2 (0, T; V^ +d ). (3.24) 

Remark 3.10. Let 49 + A9 X + 29 2 > n + 2, 29 + 29 1 > 1 - k, 9 + 29 2 > 1, 30 + 4B X > 1, 

9 + 29 x > t, and 39 + 29 1 + 29 2 > 2 - t, for some k,£e {0, 1}. Let 

13 G (^-2(9 1 +9 2 )-9,39+29 1 -^)n[^-9 1 -9 2 ,mm{29+9 2 -l,29-9 2 +29 1 -k}]. 
Then the trilinear forms b\ and 63 satisfy the hypotheses of the above theorem. 
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Remark 3.11. Let 46 + 46 1 + 26 2 > n + 2, 6 + 26 2 > 0, and 6 + 26 1 > 1. Let 

(3 G (^-2(9 1 +9 2 )-9,39+29 1 -^)f][l-9 1 -9 2 ,mm{29 + 9 2 ,29 - 9 2 + 29 1 -l}}. 

Then the trilinear form b 2 satisfies the hypotheses of the above theorem. 

Remark 3.12. Let 46 + 46 x + 26 2 > n + 2, 6 + 26 2 > 1, and 6 + 26 1 > 1. Let 

(3 G (^-2(6^6^-6, 36+26 1 -^)n[l-6 1 -6 2 ,mm{26+6 2 -l,26-6 2 + 26 1 -l}]. 

Then the trilinear forms 64 and 65 satisfy the hypotheses of the above theorem. 

Proof of Theorem 3.9. By Theorems 3. l.b) and 3. 5. a), there is s > depending on 
||w(0) such that u G L°°(0, s; V 13 ) fl L 2 (0, s; V^ +e ). With / = [0,s),wehave 

(u,A 2f3 u) + (Au,A 2/3 u) +b{u,u,A 2(3 u) = </,A 2/3 w), a.e. in /. (3.25) 

By employing the boundedness of b and the coercivity of A, we infer 

and using Gronwall's inequality, we conclude 

IK*)II?< / \\f\\l-e +K0)||Jexp f C\\u\\ 2 e _ 62 , a.e. in I, (3.26) 
Jo Jo 

where the integral in the exponent is uniformly bounded since u E L 2 (0,T; V 9 ~° 2 ). 
Therefore we have u e L°°(0,T; V p ), which transfers to u G L 2 (0,T; V d+l3 ) by the 
coercivity of A. □ 

Again for clarity, the corresponding conditions and results of Theorem 3.9 above are 
listed in the Table 5 below for the special case models listed in Table 1 . For the NS-a-like 
case in the table, the allowed values for (3 are (3 < 26 — 6 2 — 1 with /3 < 39 — |, provided 
that 6 > \ and 46 + 26 2 > 5. 

TABLE 5. Regularity results for some special cases of the model (1.1). The table 
gives values of (3 for our recovered local and global regularity results for the standard 
models, where u £ V 13 and where u € L^iV 13 ) n L 2 (F' 3+e ). (In the NS-a-like case, 
see the text for the allowed values of f3.) Again, the result for each NSE model has a 
corresponding MHD analogue. 



Model 


NSE 


Leray-a 


ML-q 


SBM 


NSV 


NS-a 


NS-a-like 


Regularity 




(3 < 1 


P<1 


(3 < 2 




(3 < 


(3 



4. Singular perturbations 

In this section, we will consider the situation where the operators A and B in the 
general three-parameter family of regularized models represented by problem (2.1) have 
values from a convergent (in a certain sense) sequence, and study the limiting behavior 
of the corresponding sequence of solutions. As special cases we have inviscid limits in 
viscous equations and limits in the a-models. 
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4. 1 . Perturbations to the linear part. Consider the problem 

ii + Au + B(u,u) = f, (4.1) 

and its perturbation 

in + A iUi + B(u h Ui) = f, (ieN), (4.2) 

where A, B, and iV (that will appear below) satisfy the assumptions stated in Section 2, 
and for i G N, : V s — > V s ~ 2e is a bounded linear operator satisfying 

\\Aiv\t £ _ e2 + \\v\\ 2 e _ e2 < (A lV ,Nv) + v G V^ 2 . (4.3) 

Assuming that both problems (4.1) and (4.2) have the same initial condition uq, and that 
Ai — > A in some topology, we are concerned with the behavior of Ui as i — ► oo. We will 
also assume that e > 9. 

Theorem 4.1. Assume the above setting, and in addition let the following conditions 
hold. 

i) b :V ai x V° 2 x 1/T->IR is bounded for some a s G [-0 2 , - 2 ]» j = 1, 2, and 
7 G [e + 6 2l (X)) n (0 2 ,oo); 

ii) &(t>, t>, iVi;) = for any v G V; 

iii) 6 : V ai x V^ 2 x V^— >R is bounded for some Oj < 9 — 9% j — 1,2, an J 7 > 7; 

iv) M G V^ 2 , an J / G L 2 (0, T; V^- 02 ), T > 0; 

v) A; converge weakly to A as i — > 00. 

TTien, fnere exists a solution u G L°°(0, T; V^ 02 ) n L 2 (0, T; V e ~° 2 ) to (4.1) such that up 
to a subsequence, 

Ui — ► w weak-star in L°°(0, T; l^ 6 * 2 ), 

— ► m wea% in L 2 (0, T; V 6 " 02 ), (4.4) 
it, — ► u strongly in L 2 (0, T; V s ) for any s < 9 — 2 , 

as z — ► 00. 

Proof. Firstly, from Theorem 3.1, we know that for z G N there exists a solution «j G 
L°°(0,T; K -6 * 2 ) n L 2 (0,T; \Z e " e2 ) to (4.2). Duality pairing (4.2) with N Ui and using 
elementary inequalities, we have 

j t (u^Nui) + 2 (A lUi ,N Ui } = 2 (f,Nui) < e^WfW 2 _ e _ 02 + e\\ Ui \\ 2 _ e2 . (4.5) 
Choosing e > small enough, then using (4.3), by Gronwall's inequality we have 

lk(*)H 2 -* 2 <e Ct . (4.6) 
Moreover, integrating (4.5), and taking into account (4.6) and (4.3), we infer 

1 1 1 1 ^ 2 (o,t; V— =— e a ) + ll' u illia(o,t;V' s - fl 2) — V'(O) ^ ^ [O7 00 )) (4.7) 

where ip : [0, 00)— >(0, 00) is a continuous function. For any fixed T > 0, this gives 
Ui G L°°(0, T; V- 02 ) n L 2 (0, T; V 9 "^) with uniformly (in i) bounded norms. On the 
other hand, we have 

||wi||- 7 < ||^i«i||- 7 + 11^(^,^)11^+ H/ll-7- (4.8) 

By estimating the second term in the right hand side as in the proof of Theorem 3.1, and 
taking into account (4.7), we conclude that u, is uniformly bounded in L 2 (0,T; V" 1 ). 
By employing Theorem A. 2, and passing to a subsequence, we infer the existence of u 
satisfying (4.4). Now taking into account the weak convergence of A4 to A, the rest of 
the proof proceeds similarly to that of Theorem 3.1. □ 
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For example, setting e = 1, with 9 = and 9 2 — 1, and checking all the requirements 
z) — f ) of Theorem 4.1, the viscous solutions to the SBM converge to the inviscid solution 
as the viscosity tends to zero. Recall that the global existence of weak solution to the 
inviscid SBM (first established in [5]) is also established in Theorem 3.1. Similarly, 
setting e = 0, with 9 = and 02 = 1, the viscous solutions to the Leray-a model 
converge to the inviscid solution as the viscosity tends to zero. This result gives another 
proof of the global existence of a weak solution for the inviscid Leray-a model. 

On the other hand, the convergence of viscous solutions of ML- a and NS-a to its 
inviscid solutions, respectively, are not covered here since both models fail condition i) 
of Theorem 4. 1 . Notice that the global existence of weak solution to the inviscid ML-a 
and NS-a are not established in Theorem 3.1. Besides the inviscid SBM, there are no 
global well-posedness results reported previously in the literature for the other inviscid 
a-models. 

4.2. Perturbations involving the nonlinear part. For i G N, let A { : V s — > V s ~ 2e and 

iVj : V s — > V s+2£2 be bounded linear operators, satisfying 

IMliU < {AiNZ\v) + \\v\\l, v G V e+e \ (4.9) 

and 

\MlZ{Nr\v), veV e \ (4.10) 

where we also assumed that iV* is invertible. In this subsection, we continue with pertur- 
bations of (4.1) of the form 

u i + A i u i + B i (u i ,u i ) = f, (ieN), (4.11) 

where Bi is some bilinear map. Again assuming that both problems (4.1) and (4.11) have 
the same initial condition u , and that Ai — > A and Bi — > B in some topology, we are 
concerned with the behavior of Ui as i — > oo. For reference, define the trilinear form 

bi(u,v,w) = (Bi(u,v),w). 

Theorem 4.2. Assume the above setting, and in addition let the following conditions 
hold. 

i) bi : V c x V c x V 7 ^IR is uniformly bounded for some a G [—9 2 , 9 + 9 2 — 2e 2 ), 
and 7 G [9 + e 2 , oo) fi (e 2 , oo); 

ii) bi(v, v, N{u) = Ofor any v G V; 

iii) M e V- Q \ and f G L 2 (0, T; V~ e ~° 2 ), T > 0; 

iv) Ai : V ~ e ' 2 — > V"~ 7 uniformly bounded and converges weakly to A; 

v) Nr 1 : V 3+202 -> y s + 2e 2- 2 ^ ^ uniformly bounded; 

vi) Nf 1 N : V /6,_6 ' 2 — >■ y e ' +6 ' 2 ~ 2e2 converges strongly to the identity map; 

vii) For a/ry t> G l/ 6 * -6 * 2 , 5j(t), t>) converges weakly to B(v, v). 

Then, there exists a solution u G L°°(0, T; V~ 02 ) n L 2 (0, T; y 9 - 02 ) to (4.1) such that up 
to a subsequence, yi = N~ x NiUi satisfies 

Ui — >• u weak-star in L°°(0, T; V~ e ' 2 ), 

Vi^u weakly in L 2 (0,T; V 9 - 02 ), (4.12) 
?/, — > u strongly in L 2 (0, T; V s ) for any s < 9 — 9 2 , 



as i — > oo. 
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Proof. Firstly, by Theorem 3.1, we know that for i £ N there exists a solution m £ 
L°°(0,T; V~ £2 ) n L 2 (0,T; V £ ~ E2 ) to (4.11). Pairing (4.11) with v t : = A^it* and using 
elementary inequalities, we have 

±{NT\,v i ) + 2(A i Nr 1 v i ,v i ) = 2(f,v i ) < e'^ff^ + e|K||^. (4.13) 

Choosing e < small enough, then using (4.9), by Gronwall's inequality and (4.10) we 
have 

\\vi(t)\\i<e Ct . (4.14) 

Moreover, integrating (4.13), and taking into account (4.14) and (4.9), we infer that for 
any fixed T > 0, v t £ L°°(0,T; V e2 ) n L 2 (0,T; r e + e2 ) with uniformly (in i) bounded 
norms. On the other hand, we have 



By estimating the right hand side as in the proof of Theorem 3.1, we conclude that it, is 
uniformly bounded in L 2 (0, T; V' 1 ), thus Vi = iVjWj is uniformly bounded in the same 
space. By employing Theorem A. 2, and passing to a subsequence, we infer the existence 

of v £ L°°(0,T; V 02 ) n L 2 (0,T; V e+e2 ) satisfying 

Vi -> u weak-star in L°°(0, T; \/ 92 ), 

«i -> u weakly in L 2 (0, T; (4.16) 
Vi —> v strongly in L 2 (0, T; \^ s ) for any s < 9 + 9 2 , 

as i — > oo. Define it = iV _1 i; and = N~ l Vi, and note that these satisfy (4.12). 

Now we will show that this limit u indeed satisfies the equation (4.1). Let w £ 
C°°(0, T; V) be an arbitrary function with w(0) = w(T) = 0. We have 

,T pT 

( Ui {t),w{t))dt+ / (AiUi{t),w{t))dt+ / 6<(tii(t), tw(t))tft 
Jo Jo 

(f(t),w(t))dt. 

We claim that each term in the above equation converges to the corresponding term in 

(u{t),w(t))dt + / (AM(t),w(t))rft + / b(u{t),u{t),w{t))dt 
Jo Jo 

T 

(f(t),w(t))dt. 



For the first term, we have 

Ui - u = N^Nyt -u = Nr l N{ yi - u) + (A^iV - I)u, 

and taking into account that Nf 1 : V s+2d2 — > V^ s+292_2e2 is uniformly bounded, and 
that N~ X N : V^ 6 *" 6 * 2 — ► \Z e+e ' 2 ~ 2e2 converges to the identity map in the strong operator 
topology, we infer Ui — > it in L 2 (0, T; \z s + 2e 2-2£ 2 ') f or an y s < 9 — 9 2 . For the second 
term, writing 

Ajitj — v4u = Ai{ui — u) + (Aj — A) u, 
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and taking into account the uniform boundedness of A h and the weak convergence A { — ► 
A, prove the claim. Finally, for the third term, we have 

Bi{ui,Ui) — B(u,u) = Bi(ui,Ui — u) + Biiui — u,u) + Bi(u,u) — B(u,u), 

and using the uniform boundedness of B { and the convergence of Bi to B, we complete 
the proof. □ 

For example, setting e = e 2 = 1, with 9 = 1 and 9 2 = 0, and checking all the 
requirements i) — vii) of Theorem 4.2, the weak solutions to the NS-a model converge 
to a weak solution of the NSE as the parameter a — > 0. This result was previously 
reported in [20] . 

5. Global attractors 

In this section we establish the existence of a global attractor for the general three- 
parameter family of regularized models, and give general requirements for estimating its 
dimension. The dimension of the global attractor gives us some measure of the level of 
complexity of the dynamics of a given flow. 

5.1. Existence of a global attractor. The following theorem establishes the existence 
of an absorbing ball in V~ 02 . Moreover, with additional conditions, it shows not only 
the existence of an absorbing ball in a higher smoothness space V@, but also that any 
solution with initial condition in V~ 02 acquires additional smoothness in an infinitesimal 
time, in particular implying that the absorbing ball in V~ 02 is compact. 

Theorem 5.1. a) Let u G L^ c (0, oo; V' 02 ) n Lf oc (0, oo; V - 02 ) be a solution to (2.11) 
with u(0) G V~ 02 . In addition, let the following conditions hold. 

(i) (Av, Nv) > c\\v\\g_g 2 for any v G V ~ 02 , with a constant c > 0. 

(ii) sup ||/||^2( t t+T . v -e-o 2) < K, where T > and K > are constants; 

Then for some constant k > and for any T' > 0, we have 

\Ht)\te 2 + \\u\\l 2{tit+T ,. v e-e 2) < e- kt \\u(0)\\-e 2 + K, t > 0, (5.1) 

where the implicit constant may depend on T' . 

b) In addition to the above hypotheses in a), for some j3 G [—6*2, 6—62] let the following 
conditions be satisfied. 

(iii) b : V 13 x V 13 x V -^R is bounded; 

(iv) (Av, (/ - Afv) > c\\v\\l +e for v G V p+e ; 

(v) sup II / 1|^2 (tjt+T!V r/>-») < K; 

(vi) ii G L^ oc (0, 00; V 13 ^ ). 
Then for any t > we have 

IK*)II? < (e- kt \H0)\\-e 2 +K)eMe- kt \H0)\\-e 2 +K), t > t , (5.2) 
where the implicit constant may depend on to- 

Remark 5.2. The trilinear forms b\ and 63 satisfy condition (iii) provided 29 > — 
26 1 - 6 2 , 29 + 2#i > 1 - k, 26 2 + 9 > 1 and 29 x + 9> k, for some k G {0, 1}. 

The trilinear formb2 satisfies condition (iii) provided 29 > r} ^ L —29\—92, 29+29\ > 1, 
29 2 + 9 > and 29 1 + B > 1. 

The trilinear forms 64 and 65 satisfy condition (iii) provided 29 > t± y l — 2B\ — 9 2 , 
29 + 29 1 > 1, 29 2 + 9 > 1 and 29 1 + 9> 1. 
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Remark 5.3. All the special cases listed in Table 1 ( except NSE) satisfy condition ( Hi). In 
the case ofNS-a-like model, condition ( Hi) is satisfied provided 9 > 1 and 9 2 > r ^ — 2. 

Proof of Theorem 5.1. We have 

j t (u,Nu) + 2(Au,Nu) = 2(f,Nu) < e^WfW l e -e 2 + 4N\\-e 2 ;e 2 \W\\ 2 e-e 2 , 
for any e > 0. By using (i) and by choosing e > sufficiently small, we get 

j t \H 2 -e 2 +c\HU 2 < 11/11 V<v (5-3) 
and since V ~° 2 V~° 2 we have 

j t M-e 2 + khf_ d2 < H/ll V, 2 , (5-4) 
with some constant k > 0. This gives 

IK*)il% 2 <e- fct ||«(0)||% 2 + /'e-^||/(T)||VfadT 

<e- fct |K0)||% 2 +K, 
and by integrating (5.3) and using (5.5), we have 

t+T Wlfc< \Ht)\te 2 + 1^ \\f\\-e-e 2 <e- kt \H0)\te 2 +K, (5.6) 
proving a). 

Now we shall prove b). As in the proof of Theorem 3.9, we get u E L 2 (0, T; V^ +e ). 
Taking w = A 2,3 u in (2.11), and using (iv) and the boundedness of b, we have 

d 

dt 

with some constant k' > 0, implying 

ll^)ll?<e- fc(t - s) exp(|| M ||i 2(M; ^ 2) )|| M (,)||J 



«III + ^NI|<||/||?_ + NIL 02 NI? ) (5.7) 



+ y s t ^ (t - T) exp(|| W ||i 2(r) ^ 2) )||/(r)||^r. 
Integrating this over s E [t — to, t] we have 

*o||w(t)||| < exp(||u||^ 2(t _ tot . ye _ 92) ) / ||u(a)||£ds 



(5.8) 



«t 

+ toexp(||M||^ 2(t _ tot . v9 _ e2) ) / \\f(T)\\l_ e dT 

J t — to 

< (e- fc 1«(0)||% 2 + + K) ex P (e- kt \\u(0)\te 2 + K), 
where we have used (5.6) and (v). This completes the proof. □ 

For example, in the case of ML-a model, conditions i) — vi) of Theorem 5.1 are 
satisfied with (5 — 0. 

In this next corollary, we combine the results in Theorem 3. 1 , Theorem 3.5 and Theo- 
rem 5.1 to show the existence of a global attractor. 

Corollary 5.4. Let the following conditions hold. 

i) b : V U1 x V 6 - 02 x V a2 ^R is bounded for some a x < 6 - 6 2 and a 2 < 6 + 6 2 
with o\ + cr 2 < 9; 
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ii) b(w, v, Nv) = Ofor any v,w E V e 02 ; 

iii) b : V* 1 x V° 2 x V^-tR is bounded for some a { < 9 - 9 2 , % = 1, % and 7 G R 

//i addition, assume that the hypotheses (i) and (iii-v) of Theorem 5.1 are satisfied. Then, 
there exists a compact attractor A <e V~ e ' 2 for the equation (2.11) which attracts the 
bounded sets ofV~ e ' 2 . Moreover, A is connected and it is the maximal bounded attractor 
in V- 02 . 

Proof. We recall that by Theorem 3.1, there exists a solution u E L^ c (0, 00; V~ 92 ) fl 
Lf oc (0 } 00; V e ~ d2 ) to (2.1 1) with any given initial data u{0) E V^ 2 . By Theorem 3.5 this 
solution is unique and depends continuously on the initial data, so we have a continuous 
semigroup S(t) : V' 02 -> V' 6 ' 2 , t > 0. Now, by Theorem 5.1 there is a ball B in V ' 02 
which is absorbing in V^° 2 , meaning that for any bounded set U C V~ 02 there exists t\ 
such that S(t)U C B for all t>t\. Therefore for any bounded set U C V~ 02 there exists 
to such that U t >t S(t)U is relatively compact in V~ d ' 2 . Finally, applying [62, Theorem 
1. 1.1.] we have that the set A = n s > U t > s S(t)B is a compact attractor for S, and the rest 
of the result is immediate. □ 

All the special cases listed in Table 1 (except NSE) satisfy the conditions of Corol- 
lary 5.4. Again, in the case of NS-a-like model, the conditions of the corollary are 
satisfied provided 9 > 1 and 9 2 > ^ - 2. 

5.2. Estimates on the dimension of the global attractor. Next we give a result which 
can be used to develop estimates on the dimension of the global attractor. To obtain 
bounds on the dimension of global attractors, we require conditions that will guarantee 
that any m-dimensional volume element in the phase space shrinks as the flow evolves. 
The general notion is that if this is the case, then the attractor can have no m-dimensional 
subsets and hence its dimension must be less than or equal to m. If one can find such 
an m < 00, then we say that the asymptotic dynamics is determined by finite number of 
degrees of freedom. 

Theorem 5.5. Let 9 > 0. Let the equation (2.11) admit the semigroup S(t) : V~ d2 — > 
V' 02 , t>0, and let X C V e ~° 2 be a bounded set such that S(t)X = X for t > 0. Let 
the following conditions hold: 

i) b : V U1 x V e ~ 62 x V**->R is bounded for some a x < 9 - 9 2 and a 2 < 9 + 9 2 
with <Ti + <j 2 < 29; 

ii) b(w, v, Nv) = Ofor any v,w E V e ~ 02 ; 

iii) For some [3 > —9 2 and p E [1, 00], 

1 f* 

e := suplimsup- / ^(r)^!!^^ < 00; (5.9) 

u £X t^oo t Jq 

iv) For some m E N, a E [0, 1), q E [l,p], and C > 0, and for any collection 
{4>i E V 9 - 02 }^ satisfying (<^, N<p k ) = 5 ik , i, k = 1, . . . , m, 

m m 

^bifaS^Nfa) < aJ2(^^i)+C\\S(t)u \\% u EX,t>0; (5.10) 

i=l i=l 

v) For any collection {fa} as above, 

{l-a)J2(A<f> i ,N<f> i )>Ce^, (5.11) 

i=l 

Then we have dn{X) < m. 
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Proof. Given u G X, the linearization of (2.11) around the solution u(t) = S(t)u , 
t > 0, is 



U + AU + B(u,U) + B(U,u) = 0, U(0) = U eV 



-02 



There exists a unique solution to the above equation, and we will denote U (t) = L(t, uq)Uo. 
One can show that for any fixed t > 0, L(t, •) : V~ 02 — > \/~ 02 is uniformly bounded on 
X, i.e., 

sup ||L(t,M O )||__ e2 ._ 02 < oo. 
Moreover one can prove that for any fixed t > 0, 

sup ii^-giy-^.o^-oii-^o aS£ ^ , 

where _D e = {(£, rj) : £, r? G X, ||£ — ?7||_e 2 < s}. Introducing the notation 

T(t,u )U = AU + B(u,U) + B(U,u), 
where u(t) = S(t)u is understood, we have 

m mm 
i=l 2=1 i=l 

q 



i=l 

implying that 

^ m m 

inf lim inf - V (T(t, u Q )4> h Nfa) > (1 - a) V] (Afc, Nfa) 



i=l i=l 



/ 1 \ * 

— C sup I limsup -IMI/3 



> 0. 

Now we apply [62, Theorem V.3.3] (see also pp. 291 therein) to complete the proof. 



□ 



Remark 5.6. Theorem 5.5 can be used to recover estimates on the dimension of the 
global attractor for the generalized model, through the application of techniques previ- 
ously used in the literature for the special cases listed in Table 1; this is a somewhat long 
calculation that we do not include here. 



6. Determining operators 

The notion of determining modes for the Navier-Stokes and MHD equations was 
first introduced in [22] as an attempt to identify and estimate the number of degrees 
of freedom in turbulent flows (cf. [15] for a thorough discussion of the role of deter- 
mining sets in turbulence theory). This concept later led to the notion of Inertial Mani- 
folds [23]. An estimate of the number of determining modes was given in [21, 41]; the 
concepts of determining nodes and determining volumes were introduced and estimated 
in [24,25,26, 39,41]. See also [40, 13, 14]. In [37, 36], a more general concept known as 
a determining operator was introduced, and the special case of determining Junctionals 
was explicitly given. 
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Following [37, 36, 35], we now define more precisely the concepts of determining 
operators and determining junctionals for weak solutions of (2.1 1). In the following two 
definitions, we consider an operator R m : V ~ 02 — > H m , where H m C H a is a finite 
dimensional subspace with some a < —9 2 . 

Definition 6.1. Let f,g& L 2 (0, oo; V~ ~ 02 ) be any two forcing functions satisfying 

hm ||/(t)-0(t)||_ 9 _fc = O, (6.1) 

t— »oo 

and let u, v G £ 2 (0, oo; \Z e - e2 ) £g corresponding solutions to (2.11). Then R m is called 
an asymptotic determining operator for weak solutions of (2.11) if 

lim \\R m [u(t)-v(t))\\-e a = 0, (6.2) 

implies that 

hm \\u(t) = 0. (6.3) 

Definition 6.2. Wfrft X C V*~ fe , letu(t),v(t) E K,t eR,be solutions to (2.11). Then 
R m is called a determining operator on the set K for weak solutions of (2.11) if 

R m u(t) = R m v(t), t e R, (6.4) 

implies that u = v. 

Given a basis i for the finite-dimensional space H m , and a set of bounded linear 

functionals {h}^ C V 02 ~ , we can construct the operator 

m 

RmU = y^Ji{u)<t>i- (6.5) 

i=l 

The assumption (6.2) is then implied by: 

lim \h(u(t) -v(t))\ =0, i = l,...,m (6.6) 

t— >tX3 

so that we can ask equivalently whether the set {h}^ forms a set of determining func- 
tionals. The analysis of whether R m or {h}^ are determining can be reduced to an 
analysis of the approximation properties of R m . Note that in this construction, the ba- 
sis {(fii}^ need not span a subspace of the solution space V~ 02 or even of H~ 02 , so 
that the functions 0j need not be divergence-free or be in H~ 02 , for example. Note that 
Definitions 6.1 and 6.2 encompasses each of the notions of determining modes, nodes, 
volumes, and functionals, by making particular choices for the sets of functions {<pi} r [L l 
and {/, };";. 

Here we extend the results of [37, 36, 35] to the generalized Navier-Stokes model 
(2.1). In particular, we will show that if {H rn C H a : m E N} is a family of finite 
dimensional subspaces, and if a family of operators R m : V 9 ~ 02 — > H rn , m G N, satisfies 
an approximation inequality of the form 

\\u-R m u\\ a < £(m)|H|0_0 2 , (6.7) 

for a function £ : (0, oo) — > (0, oo) with lim m ^oo £(m) = 0, then the operator R m is a 
determining operator in the sense of Definitions 6. 1 and 6.2, provided m is large enough. 

If H m both contains all polynomials of degree less than \6 — 62] , and is spanned by 
compactly supported functions such that the diameter of the supports is uniformly pro- 
portional to mr 1 /", we typically have £(m) ~ m -(9-92-a)/n^ p rov ided that R m realizes 
a near-best approximation of any u G V 9 ^ 02 from the subspaces H m in the if Q -norm. In 
particular, standard finite element and wavelet subspaces of sufficiently high polynomial 
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degree satisfy these conditions. Then R m may be chosen to be interpolation or quasi- 
interpolation operators, cf.[35, 37]. For example, the piecewise constants with local av- 
eraging and piecewise linears with e.g. the Scott-Zhang quasi-interpolators as in [35, 37] 
correspond to the determining volumes and determining nodes, respectively. 

Determining modes can be understood as follows. Let A : V e ~ 92 — > V a be a bound- 
edly invertible operator with the eigenvalues < Ai < A2 < • • ., such that \j — > 00 
as j — > 00. Let {wj} C V 9 ~ 92 be the corresponding set of eigenfunctions, orthonormal 
in V a . Then, with H m : = span{w 1; . . . , w m } and R m the ^"-orthogonal projector onto 
H m , it is easy to see that \\u — R m u\\ a < A" 1 ^!^-^ for any u E V 9 ~ 9 ' 2 , meaning that 
£(m) ~ A" 1 in this case. In particular, when A is a power of the Stokes operator, i.e, 
when A = (-PA)^ 2 "")/ 2 , we have A m ~ m (*-fc-«)A», and SQ _ m - {e ~e 2 -a)/r^ 
which coincides with the behavior of £ for the case in the previous paragraph. 

Bounds on the number of determining degrees of freedom are usually phrased in terms 
of a generalized Grashof number, which can be defined in the current context as 

G = Umsup||/(*)|| v - 9 - 9a . (6.8) 

t— >oo 

The definition in (6.8) generalizes the definition of Grashof number in the literature. For 
example, if the forcing term in NSE is given by / with dimensions mass x length x 
time -2 , then the nondimensional forcing / that appears in (2.1) can be defined in terms 
of /by 

/ = — 2 l (6-9) 

pv A 

where p is the density of dimensions mass per unit n-volume, v is the kinematic viscosity 
with dimensions length 2 x time -1 and L is the system size. In this case one can see that 
given a time independent forcing for NSE, the Grashof number is defined as 

T 2-1l/2 

g = — Hi/ii-i. ( 6 - 10 ) 

pv A 

It is known that if G is small enough, then the NSE possess a unique, globally stable, 
steady state solution [61]. As the Grashof number increases, the steady state goes through 
a sequence of bifurcations leading to a more complex dynamics of the flow. Hence, it is 
natural to use the Grashof number G to estimate the number of degrees of freedom of the 
solutions of the NSE as well as other turbulence models. 

6.1. Dissipative systems. In this subsection, we consider equations with 9 > 0. Note 
that Theorem 5.1 provides with examples where the conditions iii) of the following the- 
orem is satisfied (with (3 = 9 — 9 2 and p — 2). 

Theorem 6.3. (a) Let 9 > 0, and let u, v e L°°(0, 00; V' 02 ) n L 2 (0, 00; V - 02 ) be two 
weak solutions of (2.1) with the forcing functions f,gE L 2 (0, 00; V~ e 2 ), respectively. 
Let R m : V 6 ~ e ' 2 — > H m C H~ 92 , m G N, be a family of operators satisfying the 
approximation property (6.7) with a = —9 2 . In addition, with (3 > 9 — 9 2 , let the 
following conditions be fulfilled. 

i) b-.V^xV^x V a2 ^R is bounded for some o x < 9 - 9 2 and o 2 < 9 + 9 2 with 
ai + o 2 < 9; 

ii) b(w, z, Nz) = Ofor all w E V ai and z E V p ; 



1 



t+T 



iii) e := inf limsup — / ||w(r)||jgdr < 00 with p 

iv) lim \\f(t)-g(t)U-e 2 = 0; 

t^oo 
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v) lim ||.R m [w(£) — v(t)]||_ e2 = for some m satisfying 



t— too 

-2 

ca P*c a 
Then we have 

lim \\u(t) -v(t)\\-e 2 = 0. 

(b) Assume all of the above hypotheses with the time interval [0, oo) replaced by E, 
and the conditions iii)-v) are replaced by that f = g, e :— sup igM ||w(t)||^ < oo with 
p as above, and that R m u(t) = R m v(t) for all t G K with m as above. Then we have 
u = v. 

Proof Let w = u — v. Subtracting the equations (2.1) for u and v yields 

^ + Aw + B(u, u)-B(v,v) = f- g. (6. 1 1) 

at 

Pairing this with Nw, and by using condition ii), we get 

~ (w, Nw) + (Aw, Nw) = (f - g, Nw) - b(w, u, Nw) (6.12) 

Using Young's inequality, one can estimate the right hand side of equation (6.12) as 
follows: 



j t \H\ 2 N + 2 (Aw, Nw)<±\\f- g\\l e _ 62 + S\\Nf_ e2]e2 |MlU (fU3) 

+ ll & IIIMLIMI,3|ML-2£) 2 , 



where we denote \\w\\% := (w, Nw). We estimate the last term as follows: 

IMLIMMMU^ < IMI 2 ^~ A2 IMI^ 2 A2 IMI/3, (6.14) 

where Ai = ai t and A 2 = a2 ~ d2 . Using Young's inequality, we get 

, £ ii n(Ai+A 2 )g , 1 n n (2-Ai-A 2 )pm n» 

IMUMI/3lMU-2tf 2 < -\\w\\e-e 2 + ~ Hllfc HI/?- 

C[ pt 

Note that (Ai + X 2 )q = (2 — Ai — \ 2 )p = 2. Let us now choose 5 = 2 , lN p and 

e = MA, then it follows, taking into account the coercivity of A that 



d " "° " '|2 o/^f II I|2 ^ 1 II i \\2 



\W\ 



+ c A \\w\\_ e _ 02 - 2C A \\w\\_ d2 < -||/ - g\\_ e _e 2 + — ||HI-«J W II?- 



To bound the second term on the left from below, we employ the approximation assump- 
tion on R rn , which yields 



(6.15) 



|Hft+(ft(^-^-=N5jlH,-* 



1 

5 

This is of the form 



< -Jf-g\\-e- e2 + ^r 2 \\Rmw\\ 2 _ e2 . 



d 1 2 ii 1 2 



with obvious definition of x and y. 



w\\_ d2 + x\\w\\_ ga < y, 
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Lemma A.l can now be applied. Recall both ||/ — g\\-e-e 2 ~ > and ||i? m ^||_ e2 
as t — > oo by assumptions iv) and v). So taking into account iii) we have 

rt+T r-t+T 

lim / y + (r)dr = 0, limsup / x~(r)dr < oo. 

t-* 00 Jt t^oo Jt 

It remains to verify that for some T > 0, 



liminf / x(r)dr > 0. 
t-" 50 Jt 



This means we must verify the following inequality for some T > 0: 



£(m)- 2 >^ + ^limsup^ f +T \\u\\ldr. (6.16) 



-2 ^ 

ca ' pecA t^oo Jr T 

Therefore, if 



H\-e 2 + k\H\-8 2 <0, 



eW- 2 >^ + ^, (6.17) 
c A pec A 

implying that (6.16) holds for some T > 0, then by Lemma A.l, it follows that 

lim \\u{t) - v{t)\\_e 2 = lim \\w{t)\\_g 2 = 0. 

This completes the proof of (a). 

For (b), the right hand of (6.15) vanishes, implying 

d 

with some k > 0. The Gronwall inequality gives 

\\wm-B 2 < e Hs -Hw(s)\\ 2 _ e2 , 
for any t > s, and now sending s — > — oo we get the conclusion w(t) = 0. □ 

Remark 6.4. The trilinear forms b\ and b 3 satisfy the hypothesis of a) provided y + 
20i + 2 > e±2 + 2 > \, I + 20 x > k, and + 0i > if* /or some k G {0, 1}. 

The trilinear form b 2 satisfies the hypothesis of a) provided y + 20i + 20 2 > r y 2 -, 
+ 2 > 0, | + 20i > 1, and 6> + 6»! > |. 

77?e trilinear forms 64 and 65 satisfy the hypothesis of a) provided y + 20i + 2 > ^o"^, 
^ + ^2 > §, § + 20i > 1, and + 0i > \ 

Remark 6.5. From (5.1), we have e ~ ||/||y-9-e 2 — G 2 with (3 = and p = 2. Then 
condition v) of Theorem 6.3 is equivalent to the condition £(m)~ 2 > cG 2 . Assuming that 
£(m) ~ jTj-^-^-a)/" an d putting a = -6 2 , we have ^(m) ~ m' d/n . Hence m > G n/B . 

6.2. Nondissipative systems. In this subsection, we consider non-dissipative systems, 
which are represented in our generalized model when = 0. 

Theorem 6.6. (a) Let u, v G L°°(0, 00; K -6 ' 2 ) oe ftvo solutions of (2.11) with the forcing 
functions /, p G £ 2 (0, 00; V~ 02 ), respectively, and with = and 

(Av.M;) >caH% 2 , i> G V~° 2 . 

For some a < —9 2 , let R m : V~ 02 — > i? m C -f/ a , m G N, oe a family of operators 
satisfying the approximation property (6.7). In addition, with (3 > — 2 , let the following 
conditions be fulfilled. 

i) b : V a x V p x p^ 2 — >R is bounded; 

ii) 6(w, 2;, iVz) = 0/ar all w e V a and z G V 13 ; 
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iii) e := limsup ||«(t)||/3 < oo; 

iv) lim ||/(t)-<7(t)|U = 0; 

t— *oo 

v) lim ||_R m [w(t) — f(t)]||« = for some m satisfying £(m) < -rr—n-. 
t^co - - - ||0||e 

Then we have 

lim \\u{t) - v(t)||_ fc = 0. 

Assume all of the above hypotheses with the time interval [0, oo) replaced by E, 
anJ ?/ze conditions iii)-v) are replaced by that f = g, e := sup tgR ||w(t") ||jg < oo, anj 
?/za? R m u(t) = R m v(t)for all t 6 R w?Y/z m as above. Then we have u = v. 

Proof. We start as in the proof of Theorem 6.3, but instead of (6. 1 3) we get the following. 

d „ 



^ M „ ||JV + 2 (Aw, iVw) < -||/ - g\\ 2 _ e2 + 5||iV|| 2 _ e2;e2 ||w|| 2 _ e2 + \\w\\ a \\u\\p\\w\\_e 2 . 
Let us now choose 5 = ,, , ni 2 A — , then it follows, taking into account the coercivity of A 

ll iv ll_g 2; 8 2 

that 

■^\\ w \\n + c a\\w\W < ^\\f -9\\-e 2 + IWIIMUMI/j|MI-fc- 

To bound the last term from above, we employ the approximation assumption on R m , 
which yields 

<l 'H\n + (c a - t(rn)\\b\\\\u\\p) \\w\\ 2 _ d2 



dV 



< l\\f -g\\-e 2 + ll^mHUIHIIMMMI-fe- 



(6.18) 



6"' 

This is of the form 

^||w|| 2 e 2 +a:|M| 2 e2 <V + y\\w\\-e 2 , 

and an application of Lemma A.l completes the proof of (a). Part (b) is proven following 
the same argument as in the proof of Theorem 6.3(b). □ 

Remark 6.7. Note that in the case of the NSV we need — | < a < — 1 in order to satisfy 
Theorem 6.6 i) and the condition a < — 1. Choosing a = — | we get m > G 6 . This is 
consistent with the calculations in [43]. 

Remark 6.8. The trilinear forms b% and 63 satisfy the hypothesis i) of Theorem 6.6 pro- 
vided (3 + 20! + 20 2 > P + 30 2 > 1, 20i > k, and 20 x + + 2 > 1 - k, for 
some k e {0, 1}. Also, we need that a > — 29\ — (3 — 302, a > k — 20\ — 2 , and 
a>l-k-29 x - (3- 20 2 , for some k e {0, 1}. 

Remark 6.9. The trilinear form 6 2 satisfies the hypothesis i) of Theorem 6.6 provided 
(3 + 20! + 20 2 > (3 + 30 2 > 0, 20 1 > I, and 20 1 + (3 + 2 > 1. Also, we need that 
a> n±2_ 26i _p_ 3 2) a > 1 -20^02, and a>l-20 l - (3 - 20 2 . 



Remark 6.10. The trilinear forms 64 and 65 satisfy the hypothesis i) of Theorem 6.6 
provided (3 + 20 1 + 20 2 > (3 + 30 2 > 1, 20 l > 1, and 20 1 +f3 + 2 > 1. Also, we 
need that a > ^±2 - 20 x - /3 - 30 2 , a > 1 - 20 x - 2 , and a > 1 - 20i - (3 - 20 2 . 

Remark 6.11. From (5.1) e ~ ||/||v- e 2 — G. Then condition v) of Theorem 6.3 is 
equivalent to the condition £(m) -1 < cG. Assuming that £(m) ~ m~( 9 ~ 92 ~ a ^ n , and 
putting = 0, we have £(m) — m^ 2+a ^ n . Hence m > Q~ n /( e 2+ a ) t where a < —0 2 . 
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7. Length-scale estimates in terms of the Reynolds number 

In the previous section, we established estimates on the number of degrees of free- 
dom in terms of the generalized Grashof number G, a dimensionless parameter which 
measures the relative magnitude of forcing to the viscosity v. A complementary scheme 
was introduced by Doering and Foias in [17], to recast all the estimates in terms of the 
Reynolds number Re. The Reynolds number measures the effect of nonlinearity in the 
fluid response, and in the current setting it allows us to measure the effects of modifying 
the nonlinearity. It is important to recognize that in the engineering and physics com- 
munities, the Reynolds number is used more frequently than the Grashof number, as is 
viewed as more directly physically meaningful. Therefore, in this section we will derive 
a lower bound on the Kolmogorov dissipation length- scale in terms of the Reynolds num- 
ber, which will help provide some tools for relating the Grashof number-based results of 
the previous section to analogous statements involving the Reynolds number. 

To set some notation, we briefly review the ideas discussed in [17, 29] and then apply 
the analogous procedure to our more general problem. Given the velocity field u for 
the Navier-Stokes equations taken on an n-dimensional periodic domain [0, L] n with 
divergence free condition, the Reynolds number Re of the flow is defined as 

i?e = ^, where U 2 = L~ n ||u||§, (7.1) 

where the overline denotes the long-time average 

1 /■* 

g = limsup - / g(r) dr, (7.2) 



t— +oo t Jq 

and I is characterized by the following "narrow-band type" assumption on / 

|| V7||o * Hi/Ho. (7.3) 

Recall the standard definition of Grashof number in n dimensions in terms of the "root 
mean square" of the force 

C7 = ^p, where / rms = L-/ 2 ||/|| (7.4) 

Doering and Foias showed recently in [17] that in the limit G — > oo, the solutions of the 
n-dimensional Navier-Stokes equations satisfy 

G < Re 2 + Re. (7.5) 

The above estimate gives a way to transform any estimate given in terms of G into an 
analogous estimate given in terms of Re. However, as the following example from [29, 
30] shows, this procedure does not always give sharp estimates. Consider the problem of 
bounding the Kolmogorov dissipation length-scale for the Navier-Stokes equations from 
below. The Kolmogorov dissipation length-scale in terms of the energy dissipation rate 
e is given by 

l d = (z/ 3 / £ ) 1/4 , where e = uL- n \\Vu\\ 2 < G 2 . (7.6) 
Then, using (7.5) we have the following bound 

I • l d 1 < Re. (7.7) 

In the three-dimensional case, one can obtain an estimate for the number of degrees of 
freedom in turbulent flows by dividing a typical length-scale of the flow by Id and taking 
the third power. Thus, (7.7) gives an upper bound of Re 3 to the number of degrees of 
freedom in turbulent flows which is not sharp compared to the generally accepted i?e 9//4 . 
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The authors of [17, 29] obtained the bound of Re 9 ^ by time-averaging the Leray's energy 
inequality and using (7.5). 

In comparing estimates for the Navier-Stokes equations given in terms of the Reynolds 
number (and other quantities) to similar estimates for regularized equations, there is a ba- 
sic issue of identifying precisely what the Reynolds number is for the regularized equa- 
tions. Here we will extend the approach followed in [29, 30], and identify the Navier- 
Stokes velocity field as u = Nu (or u = Mu if M is more smoothing than N), where u is 
the regularized velocity field, and define the Reynolds number and the energy dissipation 
rate e in terms of u. We have then 

||u|| = || ATlfc|| ~ ||m||_26» 2 7 or INIo = ||^ M ||o ~ IMI-201) 

and similarly, ||V£t|| ~ |M|i-26» 2 or || Vw||o ~ ||u||i_2e 1 . In other words, the approach 
of [29, 30] naturally extends to our more general setting here, giving definitions of Re 
and e that satisfy 



Re 2 ~ U 2 ~ IMH 2 #> an( i £ ~ ll u lli_2e' w ^ tn = max {^i5 $2}- 

The constants in the definitions are irrelevant as far as the asymptotic behavior of the 
bounds are concerned. In the following theorem we derive a bound on a mean square 
Sobolev norm of u in terms of a similar norm with a lower Sobolev order. This will 
make possible the corollary following the theorem, which gives length-scale estimates in 
terms of the Reynolds number for the general setting in this paper. This makes it possible 
to relate the bounds of §6 giving in terms of the Grashof number to analogous bounds 
given in terms of the Reynolds number. 

Theorem 7.1. Let u e L°°(0, 00; V' 02 ) n 1^(0, 00; V - 02 ) be a solution to (2.1), and 
let a < j3 < 9 — 9 2 . Let the following conditions hold: 

i) b : V a x V a x V 7 ^IR is bounded for some 7; 

ii) b(v, v, Nv) = Ofor v e V ~ 02 ; 

iii) (Av, Nv) > c\\v\\g_Q for any v G V ~ 02 , with a constant c > 0; 

iv) The forcing term f is independent oft, and satisfies 



<C(s) 

for any s, with constants C(s) depending on s. 



Then we have 



Proof. Pairing (2.1) with Nu, and by using ii), we get 

(Au,Nu) = (f,Nu)-~(u,Nu) 

Then integrating in time gives 

t rt 



\u\ 



~ / ||w|| a ||/||_202-a + II m IIl° c ((U-V- 9 2V ( 7 - 8 ) 



6-6 2 rZ / II "Hall J \\-2V 2 ~a ~ II U \\L°°(0 t-V- 2)- 
^0 



In terms of long time averages, from the above one can derive 

MU< \\fUe 2 -«\/Ml< \\f\\oU, 



where U 2 = \\u\\^, and we have used iv) in the last step. 
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Now we will bound ||/|| in terms of u. To this end, let us pair (2.1) with /, and write 

(f,f) = (AuJ} + b(u,uJ) + j t {uj). 
Recalling that / is independent of t, and integrating in t, we have 

*ll/llo £ f + f \HlWfh + \H L ~ m v-°2)\\f\\e„ (7.9) 

Jo Jo 

which implies 

\\f\\o<U + U 2 . 
Now plugging (7.9) into (7.8), we have 

J \M\e-o, £ f NU ■ \ U \\u\\ a + jf \\u\\l + 0(1)) + 0(1), (7.10) 

giving 



\H 2 e -e 2 <U 2 + U\ 
Finally, we write the interpolation inequality 

n n ^ ii M 1 -a n 1 1 a • , x — a 

\\u\\p < \\u\\ a \\u\\ g _ g2 , withA = - — , 

u — t>2 — a 

and calculate its long time average to establish the proof. □ 

Now we apply the above result to the situation where a = —29i and — 1 — 20\, and 
the trilinear form b is given by b\ with Q\ = 6 2 . 

Corollary 7.2. Let the conditions ii), Hi) and iv) of the previous theorem be satisfied, and 
let b = bi with 6\ = 6 2 . Then we have 



e < U 2 + U 2+e + e i , where e = \\u\\l_ 2&i , and U' 2 = \\u\\ 2 _ 2Qi 
In terms of the Kolmogorov length-scale and the Reynolds number, this is 



Remark 7.3. a) For the simplified Bardina model, our result is consistent with that 
of 130], where the authors derive 

If < Re*'*. 

b) For the hyper-viscous Navier-Stokes equation and the Navier-Stokes-Voight model, 
we have the estimates 

I, 1 < Re 9 /\ and I, 1 < Re 3 / 4 , 
respectively, that appear to be new. 

8. Summary 

In this article we considered a general three-parameter family of regularized Navier- 
Stokes and MHD models on n-dimensional smooth compact Riemannian manifolds, 
with n ^ 2; this family captures most of the specific regularized models that have been 
proposed and analyzed in the literature. Well-studied members of this family include 
the Navier-Stokes equations, the Navier-Stokes-a model, the Leray-a model, the Modi- 
fied Leray-a model, the Simplified Bardina model, the Navier-Stokes-Voight model, the 
Navier-Stokes-a-like models, and several MHD models; the general model also captures 
a number of additional models that have not been specifically identified or analyzed in 
the literature. We gave a unified analysis of this entire family of models using essentially 
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only abstract mapping properties of the principal dissipation and smoothing operators, 
and then used assumptions about the specific form of the parameterizations, leading to 
specific models, only when necessary to obtain the sharpest results. In §2, we established 
our notation and gave some basic preliminary results for the operators appearing in the 
general regularized model. In §3, we built some well-posedness results for the general 
model; this included existence results (§3.1), regularity results (§3.3), and uniqueness 
and stability results (§3.2). In §4 we established some results for singular perturbations, 
which as special cases include the inviscid limit of viscous models and the a — * limit in 
a models; this involved a separate analysis of the linear (§4.1) and nonlinear (§4.2) terms. 
In §5, we showed existence of a global attractor for the general model (§5.1), and then 
gave estimates for the dimension of the global attractor (§5.2). In §6 we established some 
results on determining operators for the two distinct subfamilies of dissipative (§6.1) and 
non-dissipative (§6.2) models. In §7, we established length-scale estimates for the gen- 
eralized model; this makes it possible to recast our estimates for the number of freedom 
of turbulent flows given in §6 in terms of the Reynolds number. 

To make the paper reasonably self-contained, in Appendix A we also included some 
supporting material on Gronwall-type inequalities (Appendix A.l), spaces involving time 
(Appendix A. 2), and Sobolev spaces (Appendix A. 3). In addition to establishing a num- 
ber of technical results for all models in this general three-parameter family, the frame- 
work we developed can recover most of the existing existence, regularity, uniqueness, 
stability, attractor existence and dimension, and determining operator results for the well- 
known specific members of this family of regularized Navier-Stokes and MHD models. 
Analyzing the more abstract generalized model allows for a simpler analysis that helps 
bring out the common structure of the various models, and also helps clarify the core 
common features of many of the new and existing results. More general MHD models 
can be analyzed using the framework with only minor modifications as outlined in the 
text. 

In ongoing work, we are extending the unified analysis presented here to establish 
partial regularity results for the three -parameter generalized model, in the spirit of [3]. 
In [44], it was found that for the hyper-dissipative model, there exists a solution for 
which the Hausdorff dimension of the space-time singular set is at most 5 — AO. We 
would like to extend this result for our generalized equation to see the interplay between 
the nonlinearity, which is controlled by two parameters Oi and 2 and the dissipative term, 
which is controlled by the parameter 9 in the model equations. 

In [31] the notion of suitable weak solutions for NSE was defined. The definition in- 
troduces two parameters: a discretization scale h and a large eddy scale e. We also plan 
to extend this unified analysis to find the interplay between the nonlinear term and dissi- 
pative term that will satisfy the proposed list of mathematical criteria when establishing 
a reasonable definition of large eddy simulation (LES) models. In [3 1], the authors men- 
tioned some technical or fundamental difficulties when establishing the convergence of 
the discrete approximations of the NS-a model to suitable weak solutions of the NSE. 
We would like to use the unified analysis to see under what conditions we can recover 
the relationship between the regularizing and discretization parameters that will allow 
the model equation to be a suitable approximation to NSE. 
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Appendix A. Some key technical tools and some supporting results 

A. 1 . A Gronwall type inequality. The following is a slight generalization of the Gronwall- 
type inequality appeared in [21] and [39]. 

Lemma A.l. Let T > be fixed, and let x, y, and z be locally integrable and real-valued 
functions on (0, oo), satisfying 

rt+T rt+T rt+T 

liminf / x(r)dr > 0, limsup / x~(r)dr < oo, lim / y + (r)dr = 0, 

Jt t->°o Jt t-xnJt 

where x~ = max{— x,0} and y + = ma,x{y, 0}. If C, is an absolutely continuous non- 
negative function on (0, oo), and (, satisfies the following differential inequality 

+ x£ <y + y£ p , a.e. on (0,oo), 
for some constant p G (0, 1], then lim^oo £(£) = 0. 

A. 2. Spaces involving time. Let us recall the following well known result. A proof can 
be found in [61]. 

Theorem A.2. Let X Y <^-> Z be reflexive Banach spaces, with X <^-> Y compact. 
Let p, q > 1 be constants, and define 

dv 

y = y(0, T- p, q- X, Z) = {v G L p (0, T; X) : — G L*(0, T;Z)}. 
Then we have 3^^(7(0, T; Z) and the embedding y ■=— > L p (0, T; Y) is compact. 

A. 3. Multiplication in Sobolev spaces. With s G R, let H s be the standard Sobolev 
space on an n-dimensional compact Riemannian manifold M. We state here a well 
known result on pointwise multiplication of functions in Sobolev spaces. 

Lemma A.3. Let s, s±, and S2 be real numbers satisfying 

Si + s 2 ^ 0, Si ^ s, and si + s 2 — s > | , 

where the strictness of the last two inequalities can be interchanged ifsE No. Then, the 
pointwise multiplication of functions extends uniquely to a continuous bilinear map 

H Sl ® H S2 H s . 

Proof. A proof is given in [65] (see also [53]) for the case s > 0, and by using a duality 
argument one can easily extend the proof to negative values of s. □ 
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